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Abstract. We study the asymptotics of representations of a fixed 
' compact Lie group. We prove that the limit behavior of a sequence 

■ of such representations can be described in terms of certain random 
I matrices; in particular operations on representations (for exam- 

I— i' pie: tensor product, restriction to a subgroup) correspond to some 

, natural operations on random matrices (respectively: sum of inde- 

Oh ' pendent random matrices, taking the corners of a random matrix). 

(~| ' Our method of proof is to treat the canonical block matrix associ- 

"j^ I ated to a representation as a random matrix with non-commutative 

^ ' entries. 

(N ■ 

> : 

. 1. Introduction 

00 ■ 

^ ! 1.1. Need for an asymptotic theory of representations. One 

I i of the main questions in representation theory of Lie groups and Lie 

^ ! algebras is to understand multiplicities, commutant spaces and struc- 

"Th ! ture of representations arising in various natural situations, such as 

'. restriction to a subgroup or tensor product of representations. 

^ ! There are many reasons to study asymptotic versions of such ques- 

> ! tions in the Hmit when the representation (and, possibly, also the Lie 

! group) tends in some sense to infinity. 

■ • Prom the viewpoint of probability theory it is natural to con- 

sider the limit theorems (such as laws of large numbers, central 
Hmit theorem, etc.) in order to study the Hmits of probabil- 
ity measures on a given set. Reducible representations of a 
given group, the subject of this article, can be alternatively 
described as probability measures on the set of irreducible 
representations. 

• Even though for nearly all problems in the representation the- 
ory there are explicit answers |FH91| , IGW98j , they are based 
on some combinatorial algorithms which are too cumbersome 
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to be tractable asymptotically. For this reason in the as- 
ymptotic theory of representations one has to look for non- 
combinatorial tools such as random matrix theory |Meh91j or 
Voiculescu's free probability theory |VDN92j . 

• In theoretical physics it is a natural question as well. Indeed, 
finite dimensional representation theory is described as a nice 
discrete object which can be scaled in some thermodynamic 
limit to model continuous phenomena |KSW96a| IKSW96b] . 

• Many important questions in the theory of operator alge- 
bras concern free group factors. One of the foundations of 
Voiculescu's free probability theory was the fact that free 
products may be approximated in some sense by random ma- 
trices. It was observed by Biane |Bia95j that also representa- 
tions may provide such finite-dimensional approximants with 
an adventage of being fully constructive and non-random. 

1.2. Asymptotics of representations of a fixed Lie group. Let G 
be a fixed compact Lie group (in the following we consider the exam- 
ple when G = U(d) is the unitary group). One of our motivations is 
to study the problem of decomposing a given concrete reducible rep- 
resentation of G into irreducible components. For example we wish to 
study the following interesting examples of reducible representations: 

(1) restriction of a given irreducible representation of G' to its 
subgroup G, where G' is a given compact Lie group (for ex- 
ample: G = U(d) and G' = U(d') with d < d'); 

(2) Kronecker tensor product of given two irreducible representa- 
tions of G. 

The irreducible representations of compact Lie groups are uniquely 
determined by their highest weights (cf Section [3] for definitions); in 
the example of G = U(d] the irreducible representations are indexed 
by sequences of integers Ai > ■ ■ ■ > A^j. For explicit answers to the 
above questions there are well-known algorithms involving combina- 
torial manipulations on the the highest weights. For example, the 
decomposition of an irreducible representation of G' = U(d + 1 ) to 
a subgroup G = U(d) is given by the Weyl's branching rule; in the 
general case the answer to this question is given by counting certain 
Littelmann paths |Lit95j . However, when we are interested in the 
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situation when the dimension of the representations tends to infinity, 
these combinatorial algorithms become very cumbersome. In particu- 
lar, the direct study of the multiplicities in the above problems seems 
to be rather difficult. 

In order to avoid such difficulties we concentrate on approximate 
asymptotic answers. More explicitly, for a given sequence [pn) of rep- 
resentations of G for which the highest weights of the irreducible com- 
ponents tend to infinity, we study the limit distribution of a (rescaled) 
highest weight of a randomly chosen irreducible component of p^. 

This situation in the asymptotic theory of representions in many 
ways resembles statistical mechanics: when the number of particles 
in a physical system grows, the complexity of its description also 
grows so that its exact solution becomes quickly intractable. How- 
ever, a more modest approach in which we are interested only in 
some macroscopic quantities may result in a surprisingly simple de- 
scription. Similarly, in the asymptotic theory of representations when 
we abandon the attempts to find exact combinatorial solutions and 
restrict ourselves to a statistical description we may expect dramatic 
simplifications. Since representations have a highly non-commutative 
structure and the combinatorial algorithms behind are so cumber- 
some therefore we can expect quite surprising results. 

1.3. The main result: Representations of Lie groups and ran- 
dom matrices. For a fixed compact Lie group G and a sequence 
(Pn] of representations we construct a random matrix whose joint 
eigenvalues distribution depends on the asymptotic behavior of the 
highest weights of Pn- The symmetry of this random matrix depends 
on the group G; for example when G = U(d) is the unitary group 
the corresponding random matrix will be a hermitian d x d random 
matrix which additionally is invariant under conjugation by unitary 
matrices. 

We prove that some operations on the representations (such as 
Kronecker tensor product, restriction to a subgroup) correspond to 
some natural operations on the appropriate random matrices (sum of 
independent random matrices, taking the corner of a random matrix, 
respectively). In this way problems concerning asymptotic properties 
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of representations are reduced to much simpler analytic problems 
concerning random matrices. 

Some results proved in this paper were already considered by Heck- 
man |Hec82j . However, as far as we could understand, his proofs are 
very different from ours. Our methods are probabilistic and lead to 
many new applications and examples, see Section [5l 

1.4. The main tool: random matrices with non-commutative 
entries. Let p be a representation of a compact Lie group G; from 
the following on we shall restrict our attention to the corresponding 
representation of the Lie algebra q (if G is connected the representa- 
tion of the Lie group is uniquely determined by the representation of 
the Lie algebra). 

The family of matrices {p(x) : x G 5} can be viewed as a fam- 
ily of non- commutative random variables; in other words p can 
be viewed as a non- commutative random vector in g*. We prove 
that asymptotically, when the representation p tends to infinity, such 
a non-commutative random vector converges in distribution (after 
some rescaling) to a classical (commutative) random vector in g* 
and hence — in many cases — can be identified with a random matrix. 

This idea is closely related to the work of Kuperberg |Kup02[ 
Kup05| who — among other results — gives a new, conceptual proof 



of the result of Johansson [JohOlj (see Theorem [20] below). Kuper- 
berg's idea is to treat elements of the Lie algebra as non-commutative 
random variables and to show that for the tensor product of repre- 
sentations a non-commutative central limit theorem can be applied. 
The results of this article can be therefore viewed as an extension 
of some of Kuperberg's results |Kup02|| from central limit theorem 
related to the tensor product of representations to other operations 
on representations such as restrictions or tensor products of a fixed 
number of representations. 

1.5. Asymptotics of representations for a series of Lie groups 
and free probability. A variation of the above problem appears 
when we replace the fixed group G by a classical series of compact Lie 
groups [Gn] and we consider a series (pn), where Pn is a representation 
of Gn] we are interested in solving the analogues of the problems ([1])- 
(l2j) from Section II. 2[ For example, we may take G^ — U(n) to be 
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the series of the unitary groups. This case was studied in detail by 
Biane |Bia95j who found a connection between asymptotics of such 



representations and Voiculescu's free probability theory |VDN92j . 



In a subsequent paper |CS07j we show that the method of random 
matrices with non-commuting entries can also be applied to this sit- 
uation and the results obtained in this way are significantly stronger 
than the ones of Biane |Bia95j . 

1.6. Asymptotics of representations of symmetric groups. It 
turns out that the ideas presented in this article can be also applied 
to the case of the symmetric groups. A canonical matrix associated to 
a representation of the symmetric group was given by Biane |Bia98j 
and it turns out that the recent results of the second-named author 



|Sni06j were proved by treating (in a very concealed way) this matrix 
as a permutationally-invaria nt ran dom matrix with non-commuting 
entries. A subsequent paper (SS07 will present the details. 



1.7. Overview of this article. In Section [2] we introduce some no- 
tations concerning non-commutative random variables. In Section [3] 
we state recall some facts about representation theory and fix some 
notation. In Section [4] we prove the main result and in Section [5] we 
investigate its new consequences. 



2. NON-COMMUTATIVE RANDOM VARIABLES 

2.1. Non-commutative probability spaces. Let (n,OJt, P) be a 
Kolmogorov probability space. We consider an algebra 

Tl>1 

of random variables with all moments finite. This algebra is equipped 
with a functional E : C°°^{D.) —> M which to a random variable asso- 
ciates its mean value. 

We consider a generalization of the above setup in which the com- 
mutative algebra £°°^(0) is replaced by any (possibly non-commu- 
tative) ^-algebra 21 with a unit and E : 21 M is any linear func- 
tional which is normalized (i.e. E(l) — 1), positive (i.e. E(x*x) > 
for all X e 21 such that x 7^ 0) and tracial (i.e. E(xy) = E(yx) for 
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all X, y G 21). The elements of 21 are called non- commutative ran- 
dom variables and the functional E is called mean value or expecta- 
tion. We also say that (21, E) is a non-commutative probability space 
IVDN92I . 

Here are two motivating examples which will be used in the follow- 
ing. 

Example 1 . For any Kolmogorov probability space the corresponding 
pair (£°°^(n),E) is a non-commutative probability space. 

Example 2. For any Kolmogorov probability space and integer d > 1 
we consider the algebra C°°-[0.] » Ma = C'^^iCl;Md) of d x d 
random matrices and we equip it with a functional 

E(x) = E(trx), 

for any random matrix x, where E on the right-hand-side denotes the 
mean value and 

1 

trx = — Trx 
d 

is the normalized trace. In this way [C°°^{Cl;Md),^) is a non- 
commutative probability space. 

Any non-commutative probability space (21, E) can be equipped 
with the corresponding norm: 

||x||l2 = a/E(xx*). 

Notice that the above definitions of non-commutative probability 
spaces and random variables do not require any analytic notions other 
than positivity. In particular, as we shall see in the remaining part of 
this section, by the distribution of a non-commutative random vari- 
able we understand the collection of its mixed moments. While this 
approach turns out to to be very useful to state and prove our results 
and it has the advantage to encompass many non-bounded random 
variables, it has a couple of drawbacks. For instance, the convergence 
of non-commutative distributions as defined in Section 12.21 does not 
coincide in the commutative case with the weak convergence of prob- 
ability measures, therefore some of our corollaries concerning conver- 
gence in distribution of classical random variables will be formulated 
and proved not in the most desirable weak topology of probability 
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measures but with respect to the moments convergence. This issue 
and the way to fix it in the cases which are of our interest (so that 
the convergence in the weak topology of probability measures in fact 
holds true) are discussed in Section 14. 3[ 

2.2. Random vectors. Let V be a finite dimensional (real) vector 
space. If V : O ^ V is a random variable valued in the space V we 
say that v is a (classical) random vector in V. We say that v has 
all moments finite if E||v||'^ < oo holds true for any exponent k > 1 . 
Notice that this definition does not depend on the choice of the norm 
II ■ II on V. For a random vector v with all moments finite we define 
its moments 

(1) = m^(v) = E y ^ ■ ■ ■ O V £ V^^. 

k times 

In the case when V = M, v becomes a usual number-valued random 
variable; furthermore V®^ = R®"^ = M and the moments = Ev'^ e 
R are just real numbers and this definition coincides with the usual 
notion of the moments of a random variable. In the following we are 
interested in the space 

{v : a ^ V such that EHvH"^ < oo for each k > 1 } 

of random vectors with all moments finite which we will view as a 
tensor product 

(2) v®/:°°-(a). 

2.3. Non-commutative random vectors. Let (21, E) be anon-com- 
mutative probability space; in analogy to ([2]) we call the elements 
of V ® 21 non- commutative random vectors in V (over a non- 
commutative probability space (21, E]^. 

Given vi = xi ai e Vi (g) 21 and vi = X2 G V2 2t we define 

vi(^V2 = (xi (g) ai](§)(x2 ® az) = {x^ X2 ® aia2) G Vi (g) V2 21 

and its linear extension on non-elementary tensors. Whenever vi = vz 
with Vi = V2 one shortens the notation as v*^^ G V®^ 2t and one 
extends it by recursion to the definition of 

^gk ^ y®k ^ 
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Observe that this definition matches the definition of the tensor prod- 
uct of compact quantum groups of Woronowicz |Wor87j provided that 



21 is a quantum group and V a representation of 21. 
The \i-th order vector moment mf (v) is defined as 

m^(v) = (Id®E)v^'^e V®\ 

If there is no ambiguity we might remove the superscript E. Note 
that in the case of Example [1] when 21 is commutative the moments, 
as defined above, coincide with the old definition ([T]) of the moments 
of a random vector. 

We define the distribution of a non-commutative random vector as 
its sequence (mf (v))i^^i_2,... of moments. Accordingly, convergence in 
distribution of non-commutative random vectors is to be understood 
as convergence of the moments. 

The above definitions can be made more explicit as follows: let 
ei, . . . , Cd be a base of the vector space V. Then a (classical) random 
vector V in V can be viewed as 



(3) ^^Y- 



where at are the (random) coordinates. Then a non-commutative 
random vector can be viewed as the sum in which at are replaced 
by non-commutative random variables. One can easily see that the 
sequence of moments 

i.1 .■■■,tk 

contains nothing else but the information about the mixed moments 
of the non-commutative coordinates ai , . . . , and the convergence 
of moments is equivalent to the convergence of the mixed moments 
of Qi, . . . , ad- 

The following result provides a necessary and sufficient condition 
for a sequence of moments to be those of a commutative vector. 

Proposition 3. A non- commutative random vector v actually arises 
from a commutative probability space iff for each value of k e 
{1,2,...} the tensor m^(v] G V^"^ is invariant under the action of 
the symmetric group. 
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Proof. The necessity is trivial. For sufficiency, if (v) are invariant 
under the action by conjugation of the symmetric group, this implies 
that the GNS representation of 21 with respect to E is its abelianized 
quotient. Since E is supposed to be faithful the proof is complete. □ 



3. Preliminaries op representation theory 

3.1. Structure of compact Lie groups. In this section we recall 
some facts about Lie groups and their algebras |BtD95i[FH91j . Let G 



be a compact Lie group and let H C G be a maximal abelian subgroup 
and let f) C g be the corresponding Lie algebras. We consider the 
adjoint action of f) on the vector space q; then there is a decomposition 

(4) fl = © 0OC 

into eigenspaces. The non-zero elements a for which the correspond- 
ing eigenspace Qoc is non-trivial are called roots. We identify f)* as a 
set of functionals on q which vanish on all root spaces Qol thus 

(5) r c q\ 

Suppose that q is equipped with a G-invariant scalar product; in 
this way = 0* and p) The above isomorphisms and the inclu- 

sion f) C g allow us to consider the inclusion f)* C g*. This inclusion 
does not depend on the choice of the G-invariant scalar product on g 
and it coincides with the inclusion ([5]) considered above. More gen- 
erally, throughout all this paper, whenever we state a result related 
to g but do not mention the invariant scalar product, it means that 
the result does not depend on the choice of the scalar product. 

By definition, the Weyl group W = [g e G : gHg"^ = H} is the 
set of the elements which normalize H. The Weyl group W is always 
finite. 

Every element x G g is conjugate to some element ij G () (i.e. there 
exists an element g G G such that Adg(x) = y, where Ad : G 
End(g) denotes the adjoint representation) and y is determined only 
up to the action of the Weyl group. Similarly, every element x G g* 
is conjugate to some element y G f)* (which is determined only up to 
the action of the Weyl group). 
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Example 4. For G = U(n) we may take H to be the group of diagonal 
unitary matrices. Then g is the set of n x n antihermitian matrices 
and is the set of diagonal matrices with imaginary entries. In this 
case W is the group of permutation matrices and hence is isomorphic 
to the symmetric group Sn- We equip q with a G-invariant scalar 
product {x,y) — Trxy* and thus we identify Q = Q* and i) = i)*. 
Clearly, every element x G = g* is conjugate by a unitary matrix 
to some diagonal matrix y e i) = I)* which is determined by the 
eigenvalues of x. The freedom of choosing the order of the diagonal 
elements of y corresponds to the action of the Weyl group S^. 

Remark 5. In random matrix theory it is more customary to work 
with the space of hermitian matrices instead of the space of antihermi- 
tian matrices; for this reason we may consider a simple isomorphism 
between these spaces given by multiplication by i. 

The above example heuristically motivates that we may think that 
for a given element x e q* the corresponding element y & i)* contains 
the information about the eigenvalues of x. 

3.2. Irreducible representations and highest weights. Let p : 
— > End(V) be an irreducible representation. The representation 
space can be decomposed into the eigenspaces of the maximal abelian 
subalgebra f): 

(6) V = V«, 

in other words H e f) acts on as multiplication by cx(h.). The 
elements a for which the corresponding eigenspace Va is non- 
trivial are called weights. 

Let us fix some element t e I) which is generic in a sense that for 
any weight non-zero oc we have a(t) 0. Now a weight a is called 
positive if a(t) > 0, otherwise it is called negative; similar conven- 
tion concerns roots as well. The highest weight of a representation 
is the weight a for which a(t) takes the maximal value. In fact, any 
irreducible representation is uniquely determined (up to equivalence) 
by its highest weight. The weight space 0oc corresponding to the high- 
est weight is always one-dimensional. Clearly, the above definitions 
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depend on the choice of t, nevertheless the notion of highest weight 
is well-defined (up to the action of the Weyl group). 

3.3. Enveloping algebra. The enveloping algebra of g, denoted by 
Q, is the free algebra generated by the elements of q quotiented by 
relations gh, — Kg — [g,h] for any g, h. G g. Usually we work over the 
complex (complexified) algebra. This algebra is naturally endowed 
with a ^-algebra structure obtained by declaring that the elements 
of g (before complexification) are anti-hermitian elements. For our 
purposes, it has the following two important features: 

• an irreducible representation (p, V) of g gives rise to an onto 
algebra homomorphism g — > End(V] and conversely; 

• therefore, irreducible representations of g are in one to one 
correspondence with minimal tracial states on g. 

We will need the following theorem, known as Poincare-Birkhoff- 
Witt theorem. 

Theorem 6. If h^, . . .hn is a basis of g as a vector space, then 

is a basis of g as a vector space. In particular, there is a filtration 
on g defined as follows: the degree ofipegis the smallest k such 
that ip is a sum of monomials of elements in g with at most k 
factors. 

3.4. Reducible representations and random highest weights. 
Let p : G ^ End(V] be a (possibly reducible) representation of G on 
a finite-dimensional vector space V. We may decompose p as a sum 
of irreducible representations: 

P = ^aPa, 

where Pa denotes the irreducible representation of G with the highest 
weight A and ua G {0, 1 ,2, . . .) denotes its multipUcity in p. 

We define a probability measure |a.p on i)* such that the probability 
of A e {)* is equal to 

Ua ■ (dimension of Pa] 
(dimension of V) 
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in other words it is proportional to the total dimension of all the 
summands of type [px] in p. In this way the probability measure |j.p 
encodes in a compact way the information about the decomposition of 
p into irreducible components. We define the random highest weight 
associated to representation p as a random variable distributed ac- 
cordingly. 

Remark 7. If we literally follow the above definition then \ip is a 
certain probability measure on the Weyl chamber. This definition has 
a disadvantage that it depends on the choice of the Weyl chamber, 
therefore sometimes it will be convenient to understand by |j,p the 
W-invariant probability measure on i)* obtained by symmetrizing the 
above measure by the action of the Weyl group W. 

3.5. Random matrix with specified eigenvalues. For an element 
A G ()* we consider a random vector in q* given by A = AdgA G g*, 
where g is a random element of G, distributed according to the Haar 
measure on G. We will say that A is a G-invariant random matrix 
with the eigenvalues A. This terminology was motivated by the 
following example. 

Example 8. If G U(n) is the group of the unitary matrices and A = 
diag(Ai , . . . , An) G [)* is a diagonal matrix then the distribution of A is 
indeed equal to the uniform measure on the manifold of antihermitian 
matrices with the eigenvalues Ai, . . . , An. 

We may extend the above definition to the case when A G [)* is 
random; in such a case we will additionally assume that g and A are 
independent. If |j. is a probability measure on f)* we may treat it as 
distribution of the random variable A; in this case we we say that A is 
a G-invariant random matrix with the distribution of eigenvalues 
given by \i. We denote by p. the corresponding distribution of A. 

Remark 9. It would be more appropriate to call A a G-invariant ran- 
dom vector in g*, nevertheless in the most interesting examples the 
Lie groups under consideration carry some canonical matrix structure 
hence the elements of g* can be indeed viewed as matrices. 

The moments of G-invariant matrices are characterized by the fol- 
lowing lemma. 
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Lemma 10. Let [i be a probability measure on f)* with all moments 
finite and let p. be the corresponding distribution of a G -invariant 
random matrix. 

For each k the moment mic = mk(p.) G (g*)*^^ is the unique 
element such that: 

(1) TTLic is invariant under the adjoint action of G, 

(2) actions of : Q^^ —> M and midix) : g®^ M coincide on 
G -invariant tensors mg®^ (where above mic(|J-) denotes the 
canonical extension of m]c{[i] E (p)*)®'^ which is possible 
thanks to the inclusion ()* C g*^. 

Proof. Point follows easily from the invariance of the Haar mea- 
sure. 

For any elements xi,X2 G g* which are conjugate to each other the 
restrictions of the maps xf^, xf ^ : g®'^ ^ R to G-invariant tensors 
coincide. In this way point (l2l) follows easily. 

For X G let 



x' = 



Adf (x) dg 



JG 

be the average over the Haar measure on G. Clearly, x' is a G- 
invariant tensor. Since 

(mk(p,))(x) = (mk(p,))(x') 

therefore the values of the functional mic(p.) are uniquely determined 
by its values on G-invariant tensors which shows the uniqueness of 

mic. □ 

Remark 1 1 . Similarly as in Remark [9l it is sometimes more conve- 
nient to understand the distribution of the eigenvalues of a random 
element of g* as a W-invariant measure on f)*. 



4. Representations and random matrices with 
non-commutative entries 

4.1. The main result. Let p : g — > End(V) be a representation of 
g. We shall view p as an element of g* ® End(V]; in other words p 
is a non-commutative random vector in g* over a non-commutative 
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probability space (End(V),tr). The sequence of its moments mi^ = 
TT^k(p) £ (0'^)®'^, or equivalently, vrik : 0®'^ ^ C is given explicitly by 

Tni<(gi (g) • • ■ (g) g^,) = tr [p(gO ■ ■ ■ p(gk)] . 

The following theorem is the main result of this article. 

Theorem 12. Let [en] be a sequence of real numbers which con- 
verges to zero. For each n let Pn '■ Q End(Vn,) be a representa- 
tion of g and let 'Kn be the corresponding random highest weight. 
Let A be a G -invariant random matrix with all moments finite. 
Then the following conditions are equivalent: 

(1) the distributions of the random variables enAn converge in 
moments to the distribution of the eigenvalues of A; 

(2) the sequence CnPn of non- commutative random matrices 
converges in distribution to A. 

Proof. Suppose that condition ([T]) holds true. It is enough to prove 
that from any subsequence (ek(n)Pk(n)) one can choose a subsequence 
(ek(i(n))Pk{i(n))) which converges in distribution to the random matrix 
A. 

Let a subsequence (ek(n)Pk{n)) be given; by LemmaflSland the com- 
pactness argument it follows that there exists a subsequence (ek{i(n))Pk{i(n))] 
which converges in distribution to some random matrix M G ® 21 
with non-commutative entries (for some non-commutative probabil- 
ity space (2t, E)). It is enough to prove that M is the G-invariant 
random matrix in q* with the same distribution of the eigenvalues as 
for A. In order to keep the notation simple instead of (ek(i(n))Pk(i(n))) 
we will write [enPn]- 

Firstly, observe that Lemma fT3] shows that 

||[enPn(Xi), enPn(X2)]||L2 = 1 1 enPn( [^1 , xj) 1 1 < 

(E||enAnf)" en||[xi,X2]|| ^ 

where we used the fact that the first factor on the right-hand side con- 
verges to some constant depending on the distribution of eigenvalues 
of A. This shows that the elements {M.(x) G 21 : x G g} commute 
hence M can be identified with a classical random variable (valued 
in q). 
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Secondly, since each of the random matrices enPn is G-invariant, 
the random matrix M. is also G-invariant. 

Thirdly, let x E g®"^ be invariant under the action of Ad®^ : G 
End(0'®^). For every irreducible representation p we have that p'^^[x) 
is a multiple of identity hence can be identified with a complex num- 
ber. The exact value of this number is equal to p'®'^(x)|y^ which can 
be estimated with the help of Lemma [HI Therefore for a (possibly 
reducible) representation Pn we have 

(mk(enpn)) M = IEA®'^(x] + (terms of degree at least 1 in e), 

where A^t^ = CnAn, where An is the random highest weight associated 
to the representation p^. Hence 

{mu{M))[x] = lim (m^lenPn)) (x) = 

lim (mk(enAn))(x) = (mk(A))(x). 

n— >oo 

The above equality and Lemma [10] show that the distribution of the 
eigenvalues of M. coincides with the distribution of eigenvalues of A 
which finishes the proof. 

Suppose that condition ([2]) holds true. In order to prove — 
similarly as in the proof of the implication ([1]) ([2]) — it is enough 
to show that from any subsequence (eic(n)Ai<{rL)) one can choose a sub- 
sequence (ek(i(rL))Ak(i(n))) which converges in moments to the distribu- 
tion of the eigenvalues of A. 

Lemma [T5l can be used to show that if k is even then the sequence 
of moments m^^ien^n) is bounded hence for every k the sequence of 
moments m^^len^n) is bounded. Again, by a compactness argument 
we can find a subsequence which converges in moments to the distri- 
bution of eigenvalues of some random matrix A'; from the implication 
([ID =^ ([2D it follows that random matrices A and A' must have equal 
moments of their entries. Prom Lemma [TBI it follows that the eigen- 
values of A' have the same moments as the eigenvalues of A which 
finishes the proof. □ 

4.2. Key lemmas. We start with the following estimate: 

Lemma 13. We equip g with a G-invariant scalar product; in this 
way we equip g = g* with the corresponding norm. 
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Let a unitary representation of a group G on a finite- dimen- 
sional Hilbert space be given and let p be the corresponding rep- 
resentation of a Lie algebra q. If p is irreducible with highest 
weight A then for any x G g 

(7) IIpMII < ||x|| ||A||, 

where the norm on the left-hand side denotes the operator norm. 

Proof. Since x e g is conjugate to some element of f), it is enough to 
prove (I7l) for X G f). For such elements the action of p(x) is diagonal 
with respect to the decomposition ([6l) and 

||p(x)|| = max|a(x]| < |A(x)|, 

a. 

where the maximum runs over the set of roots a G g* contributing to 

For a representation p : g ^ End(V) we define p^ : g®'^ End(V) 
on simple tensors by 

P^(gi ® gic) = p(gi]---p(gk) 
and extend it to the general case by linearity. 

Lemma 14. Let e be a number, Zet p : g — > End(V) be an irre- 
ducible representation of q and let A be the corresponding highest 
weight. 

Let z G g®^ be given; by a small abuse of notation we will denote 
by p^(z)|y^ : Va ^ Va the restriction of p^{z) to Va projected again 
onto Va. Since the highest- weight space Va in the decomposition 
(l6l) is one- dimensional, we will identify p^{z)\^^ with a complex 
number. Furthermore, e'^p'^(z)|y^ is a polynomial in {A(x) : x G g} 
and e, where A = eA. This polynomial can be written as 

e'^p^(z)|y^ = A®'^(z) + [terms of degree at least 1 in e). 
Proof. It is enough to prove that 

e^p(gi) ■ • • p(gic)|v, = 

A(gi) • ■ ■ A(gk) + (terms of degree at least 1 in e) if gi, . . . , g^ G J), 
(terms of degree at least 1 in e) otherwise 
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holds true for all tuples gi , . . . , gic G such that each gt belongs either 
to f) or to one of the root spaces in the decomposition ([11). We use 
induction over k: assume that the lemma holds true for all k' < k. 

Let 7r G Sk be a permutation. Thanks to the commutation relations 
xy =yx + [x, y] in the universal enveloping algebra of g we may write 
in the universal enveloping algebra: 

gi ■ • ■ gk = 9n{^) ■ ■ ■ Qnik] + (summauds with at most k — 1 factors) 

hence 

e^p(gi)- ■ ■Pi9k)\y^ = e''p(g7t(i))- ■ ■P(g7t(k))|v^ + 



e p(summands with at most k— 1 factors) 



V, 



A 



The induction hypothesis can be applied to the second summand on 
the right-hand side and it shows that it is of degree at least 1 in e, 
therefore 

e^p(gi)- ■ •p(gk)|v-, = e''p(g7t(i))- • •p(g7t(k))|v-,+ 

(terms of degree at least 1 in e). 

It follows that it is enough to consider the case when the elements 
gi , . . . , gi< are sorted in such a way that for some r, s, t > such that 
r + s + t = k the initial r elements gi, . . . , gr belong to root spaces 
corresponding to the negative roots, the next s elements gr+i , . . . , gr+s 
belong to f) and the final t elements gr+s+i , • • • , gr+s+t belong to root 
spaces corresponding to the positive roots. A direct calculation shows 
that 

e^Mg^)■■■^Q^,) if gi, . . . , g^ G f), 
otherwise 



e^P(gi) ■ ■■p[9^c] 



IVa 

and thus the inductive step follows. □ 

Lemma 15. For a given irreducible representation of q we 
denote 

MA=-}^p^x,)^ 

i 

where (xt) denotes an orthogonal basis of g (regarded as a real 
vector space). is a multiple of identity hence can be identified 
with a complex number. 
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There exists a constant C with a property that 

|A|^<2Ma + C 

for any value o/A. 

Proof. Let (Ci) be some linear basis of g (this time regarded as a 
complex vector space) and (fi) be its dual base. Then 

MA = -}^p"(eOp"(fO. 

i 

It will be convenient for us to take as (Ci) a union of two families: 
firstly from each root space Qa we select some non-zero vector (the 
corresponding dual vector fi belongs to g_cc) a-^d secondly we select 
some base of P). 

Since Ma is a multiple of identity Lemma[l4]can be used to evaluate 
it. It is easy to check that there is some element x G f) (which does 
not depend on the choice of A) with a property that 

Ma = |A|^ + A(x). 

The estimate 

2A(x) > -|A|^-|x|^ 
finishes the proof. □ 

Lemma 16. Let A be a random vector in f) which is invariant 
under the action of the Weyl group W and let A be the corre- 
sponding G -invariant random vector in g. Then each moment 
m)c(A) is a polynomial function in the moment mk(A]. 

Proof. Let us assume for simplicity that g is semisimple. The Harish- 
Chandra isomorphism (see |Kna02j for a reference) is an isomorphism 
between Z[q) (the center of the enveloping algebra g) and S(f))^ (the 
W-invariant part of the symmetric algebra S ( f) ) ) . 

The isomorphism of vector spaces g = Syin^(g) implies that 

(8) Z(g)=0[Sym^(g)]^; 

k>0 

similarly 

(9) S([))^-0[Sym^([))]^. 

k>0 
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The Harish-Chandra isomorphism provides an isomorphism between 
the right-hand sides of ([8]) and ([9l). Since it preserves the gradation 
it is the required map. 

The general case follows from the fact that the Lie algebra g can be 
written as a direct sum of its center and a semisimple Lie algebra. □ 

4.3. Weak topology and uniqueness. In view of Theorem [T2l it is 
natural to address the question whether the notion of convergence 
of moments considered there could be replaced by some more proba- 
bilistic notion of convergence. 

The notion of convergence in moments is very effective in the 
framework of bounded operators, as it can be seen for example in 
Voiculescu's free probability theory |VDN92j . Indeed, the conver- 



gence in moments fully determines the von Neumann algebra gen- 
erated by the limiting operator, and in the commutative case the 
moments of a bounded random variable determine its distribution by 
Stone-Weierstrass theorem. 

Unfortunately, if the random variables under consideration are not 
bounded, then their moments might not be finite; even if the latter 
case holds then in general the moments do not determine the dis- 
tribution of a random variable. Therefore it would be desirable to 
use some more refined description of the joint distribution of ran- 
dom variables. In the context of classical probability theory such a 
description is given by an appropriate probability measure; unfortu- 
nately, it is not clear what would be a good notion of distribution of 
unbounded operators in the non-commutative case. Some attempts 
to define weak convergence of joint distribution in the context of 
W*-probability spaces have been made in very specific examples (see 
|Mey93| and references therein). Unfortunately, it is not clear to 
us how one can adapt these definitions in our setting. A promising 
approach to this problem via Gromov-Hausdorff distance was pre- 
sented by Rieffel |Rie04j , however for the moment it is not clear if 



this method can be succesfully used for our purposes. 

To summarize the above discussion: we have no candidate for some 
kind of convergence which would replace the convergence in moments 
in point ([2l) of Theorem [12] since we deal here with a joint distribution 
of a family of non-commuting random variables. 
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Nevertheless, in point 1^ of Theorem [12] we deal with classical 
random variables therefore it makes sense to consider convergence in 
some other sense, such as weak convergence of probability measures. 
Unfortunately, in general there is no connection between convergence 
in moments and weak convergence of probability measures. In partic- 
ular, in order for convergence of measures in moments to imply their 
weak convergence we need to assume, for example, that the limit 
measure is uniquely determined by its moments. 

It seems to be hard to incorporate the weak convergence of proba- 
bility measures to condition ([T]) and preserve the equivalence of con- 
ditions and ([2]). Therefore our strategy will be to keep Theorem 
[T2] unchanged and in the study of its applications to pay attention 
to the weak convergence of probability measures (e.g. Theorem [171 
item (2) and Theorem [IH item (2)). 

5. Applications op the main theorem 

5.1. Restriction of representations and tensor product of rep- 
resentations. In this section we investigate a few remarkable conse- 
quences of Theorem [T2l 

Theorem 17. Let G C G' fee Lie groups and g C g' be the cor- 
responding Lie algebras, let (Cn) be a sequence of real numbers 
which converges to zero. Let (p^) be a sequence of representa- 
tions of g' ; by and we denote the random highest weight 
corresponding to representations and p'Jg, respectively. 

(1) Assume that CnA^^ converges in moments to the distribu- 
tion of a G' -invariant random vector A with values in g' . 
Then the sequence CnAn converges in moments towards the 
G-invariant random vector Tlg[A) where Tl^ : g' —> g is the 
orthogonal projection. 

(2) Assume that CnA^,^ converges weakly to the distribution of a 
G' -invariant random vector A with values in g' . Then the 
sequence e-nXn converges weakly towards the G-invariant 
random vector TTg(A]. 

Proof. Notice that the non-commutative random vector p'Jg is a pro- 
jection of the non-commutative random vector p'^ onto g. It follows 
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that the random matrix which is the entry wise limit of enP^^^lg is a 
projection of A to g. Theorem [12] can be applied twice: for the se- 
quence (p^) and for the sequence [Pnlg] which finishes the proof of 
part (p. 

In order to prove ([2]) it is enough to show that for every £ > and 
every subsequence eic(n)Ak{TL) we can chose a subsequence eic(i(TL))Ai<(i(rL)) 
which converges weakly to some limit distribution with a property 
that its variation distance from the distribution of ng(A) is smaller 
than £. 

Let £ > and a subsequence ek(TL)Ak(n) be fixed. For simplicity, in 
the following instead of e^ri)^^{n) we shall consider just the sequence 
enAn- We can find a sequence of representations ( p^) for which cor- 
responding rescaled random highest weights CnA^ have a common 
compact support and the total variation distance between the dis- 
tribution of GnA^ and e-J^n is smaller than e (such a sequence p'^ 
can be constructed by truncating the distribution of enA^^ to some 
sufficiently big compact set). 

By compactness argument we can select a subsequence ei{TL)A((^j 
which converges weakly (hence in moments) to some limit; let A be a 
G-invariant random vector in g* with this distribution of eigenvalues. 
The total variation distance between the distribution of A and A is 
bounded by £. 

The first part of the theorem can be applied to the subsequence of 
representations ei(TL)Pi(n)i it follows that the highest weights ei,(n)Ai,{n) 
(corresponding to the restrictions of p^n) to g) converge in moments 
to the distribution of the eigenvalues of TTg(A). Since the random 
weights have a common compact support the convergence holds also 
in the weak sense. We finish the proof by observing that the total 
variation distance between the distribution of TTglA) and TTgtA] is 
smaller than £. □ 

In the case of the inclusion of the groups U(d) C U(d') for d < d' 
the above theorem takes the following concrete form: 

Corollary 18. Let d < d' fee positive integers and let (en) be 
a sequence of real numbers which converges to zero. Let A — 
(Aij)i<ij<d' be a hermitian U.{d) -invariant random matrix and let 
Pn be a sequence of representations of U(d') with a property that 
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the distribution of of en^'n converges to the joint distribution of 
eigenvalues of A, where G Z^' is a random weight associated 
to pn. 

Then the distribution of of e-a^n converges to the joint distri- 
bution of eigenvalues of the corner (Aij)i<t j<d, where An G is 
a random weight associated to the restriction enPnIud • 

Observe that in the above result we used Remark [5] in order to 
work with hermitian random matrices. Similar concrete interpreta- 
tion in the case of unitary groups is also possible for the other results 
presented in this section. 

We leave to the Reader to investigate the other such simple inclu- 
sions: of orthogonal groups 0(d) C O(d') and of symplectic groups 
Sp(d)cSp(d'). 

Theorem 19. Let {e-n.) be a sequence of real numbers which con- 
verges to zero, let (pn'), (Pn') be two sequences of representations 
of Q and let\n^ andXn^ be the corresponding sequences of random 
highest weights. Furthermore, let An be the sequence of random 
highest weights corresponding to the tensor products 
Let A'^', A'^' be independent G-invariant random matrices in g*. 

(1) If for each i G {1,2} the sequence CnAn' converges in mo- 
ments to the distribution of eigenvalues of A'^' then CnAn 
converges in moments to the distribution of eigenvalues of 

An) + A(2). 

(2) If for each i G {1,2} the sequence enAn^ converges weakly to 
the distribution of eigenvalues o/A'^' then enAn converges 
weakly to the distribution of eigenvalues o/ A'^' + A'^'. 

Proof Let pit^' plx^' ® p^^K Then 

(10) p^^\x) = p^^\x)(^^ + ^(^Q^^\x). 

It follows that CnPn' viewed as a non-commutative random vector 
is a sum of two non-commutative random vectors: (enPn'(x)) ® 1 
and 1 ® (enPn'(x)). Theorem [T2] implies that the first summand 
converges in moments to A'^' and the second summand converges 
in moments to A*^'. The coordinates of the first vector (viewed as 
non-commutative random variables) commute with the coordinates 
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of the second vector; since we consider them with respect to a state 
(i.e. the normalized trace) which is a tensor product of the original 
states (i.e. normalized traces), it follows that their sum converges to 
the sum of the independent random matrices which finishes the proof 
of part (p. 

Part ([2]) follows from part ([1]) in a similar way as in Theorem [171 □ 

5.2. Central limit theorem. The following theorem was already 
proved by Kuperberg |Kup02| in a slightly different setting. We in- 
clude this theorem here because we believe that the proof in this 
framework (using Theorem [l2|) is new. 

Theorem 20. Let G be a Lie group and let p be its representation. 
We define c E g* given by c(x) = tr p(x). 

Let An be the random highest weight corresponding to the rep- 
resentation p*^"^. Then the sequence 

(11) _L[An-cn] 

converges (both in moments and weakly) to the distribution of 
eigenvalues of a certain centered Gaussian random matrix in q*. 
The covariance of the above Gaussian random matrix is given 

by 

A(x) A(^J)d^(A) = trp(x)p(ij) 

for any x,y e g. 

Proof. Similarly as in Eq. (fTOf l we have that 

p®^(x) = p(x) ® 1 ® ■ • • + 1 ® p(x) ® 1 ® ••• + ••• . 

regarded as a non-commutative random variable is a sum of n com- 
muting summands therefore the non-commutative central limit the- 
orem of Giri and von Waldenfels |GvW78j can be applied. It fol- 
lows that the distribution of the non-commutative random vector 
_^|^p®^_ cn] converges in moments to the distribution of a certain 
centered Gaussian random variable X which takes values in g*. We ap- 
ply Theorem [121 it follows that the distribution of the random weight 
(flTI l converges in moments to the distribution of the eigenvalues of 
X. 
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Let II ■ II be the norm on g* associated to any G-invariant scalar 
product. The distribution of X is multidimensional Gaussian there- 
fore there exists a constant d > such that the tail estimate 

(12) P(||X|| >t) < e-'^^' 

holds true for sufiicienly big values of t. For any element of its 
norm is equal to the norm of the element of [)* corresponding to 
its eigenvalues; it follows that the estimate (IT2l l remains true if the 
random variable X is replaced by its eigenvalues. The estimate ( 1121 ) 
shows therefore that the even moments of the eigenvalues distribution 
of X are dominated by the even moments of a Gaussian distribution. 

This implies that E(e'l^ll) < oo and by Corollary 2.2 and Theorem 
2.1 of |DLY02j it follows that the distribution of eigenvalues of X 
is uniquely determined by its moments. It follows by a standard 
compactness argument that the distribution of the random weights 
(flTI l converges weakly to the distribution of the eigenvalues of X. □ 

Often we have some additional information about the considered 
Lie group G which restricts the number of G-invariant Gaussian mea- 
sures on g*. In the case of G = U(d) it is convenient to consider a 
centered hermitian Gaussian random matrix g = (gij)i<i,j<d defined 
by the convariance 

Egijgki = 0, Eg^gl^ = SuSjk", 

this kind of random matrix (and the corresponding measure on the 
space of d X d hermitian matrices) is called Gaussian Unitary En- 
semble (GUE) and plays an important role in random matrix theory. 
For any v > we define GUEv as the distribution of a random matrix 

g - tr g + X 

where x is an independent centered Gaussian variable with the vari- 
ance V. 

Under the isomorphism from Remark [H GUEv becomes a measure 
on the Lie algebra u(d) and it is not very difficult to check that 
(except for degenerate cases) every U(d)-invariant Gaussian measure 
(up to dilation by some number) is of this form. In particular, we get 
the following result. 
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Corollary 21. Let p be a representation of U(d). There exist 
constants Ci,C2 with a property that if \n = (An,i , • • • , ArL,d) is a 
random weight associated to p®"^ then the joint distribution of 
the components of the vector 

-^(An-nC2j 

converges to the joint distribution of the eigenvalues of the GUEy 
random matrix. 

5.3. Toy example: representations of SO (3) and SU(2). The 
above results in the simplest non-trivial case of G = S0(3) should 
not be very surprising from a viewpoint of quantum mechanics. Each 
quantum-mechanical system in three-dimensional space can be viewed 
as a (possibly reducible) representation of SO (3) (or its universal 
cover Spin(3) = SU(2)) on some Hilbert space V. The irreducible 
components of this representation have a nice physical interpretation 
as physical states with a well-defined length |J| of the angular momen- 
tum. For simplicity, we assume that V itself is irreducible hence V is 
finite-dimensional. The information about the state of the physical 
system is encoded by a state 4> on the algebra generated by observ- 
ables. We are interested in the situation when the physical state of 
the system is SO (3)-invariant; it follows that (j) is the normalized 
trace on End(V). 

The physicist's question about the distribution of a component Jz 
of the angular momentum can be reformulated in the language of 
mathematics as question about the decomposition into irreducible 
components of the restriction V iso(2] °^ representation V to a 
subgroup SO (2) (or, more generally, restriction of the representation 
V of Spin(3) = SU(2) to its subgroup Spin(2) = U(l)), namely it is 
the uniform measure on the set of integers (or half-integers) 

(13) {-III, -IJI + 1,. ..,111-1,111} 

(for simplicity we use the system of units in which the Planck's con- 
stant ?i = 1). 

On the other side it is well-known that when the size of our sys- 
tem becomes macroscopic then quantum mechanics may be approx- 
imated by the classical mechanics, where the angular momentum 
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J = (Jx, Jy, Jz) is just a usual vector consisting of numbers. Our as- 
sumptions on irreducibility and S0(3]-invariance imply that in the 
classical limit J is a random vector with a uniform distribution on a 
sphere of a fixed length |J|. The physicists question about the distri- 
bution of a component of the angular momentum in this context 
is answered by the theorem of Archimedes, namely it is the uniform 
distribution on the interval [— |J|,|J|]. In the limit |J| — > oo, after 
appropriate rescaling, the uniform measure on the set ( 1131 ) indeed 
converges to the uniform measure on the set [— |JUJ|] hence the an- 
swer given by classical mechanics is indeed the limit of the answer 
given by quantum mechanics. 

It is also a consequence of this article: indeed we view the represen- 
tation p of the Lie algebra so(3) as an element of (so(3))* ® End(V); 
in our current case this takes a concrete form of a matrix 



(14) 



Jz -Jy 
-Jz h 
h -Jx 



(the form of this matrix depends on particular choice of the identi- 
fication of (so(3))* with certain 3x3 matrices). Theorem [T2] states 
that asymptotically the matrix (fT4f) behaves like a random 3x3 an- 
tisymmetric matrix with eigenvalues 0, |J|, — | J|. 

We leave the analysis of the central limit theorem in this case to 
the Reader. 
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